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ABSTRACT 



We present Spitzer IRAC images, along with representative 5.27 to 
38.5 jum IRS spectra of the Cassiopeia A supernova remnant. We find that 
various IRAC channels are each most sensitive to a different spectral and phys- 
ical component. Channel 1 (3.6 /tni) provides an excellent match to the radio 
synchrotron images. Where Channel 1 is strong with respect to the other IRAC 
channels, the longer-wavelength spectra show a broad continuum gently peaking 
around 26 /im, with weak or no lines. We suggest that this is due to un-enriched 
progenitor circumstellar dust behind the outer shock, heated and potentially pro- 
cessed by the photons and electrons from the shock. Where Channel 4 (8 /im) 
is bright relative to the other IRAC channels, the long-wavelength spectra show 
a strong, 2-3 /im-wide peak at 21 /im, likely due to silicates and proto-silicates. 
Strong ionic lines of [Ar II], [Ar III], [S IV] and [Ne II] also appear in these strong 
Channel 4 regions. We suggest that in these locations, the dust and ionic emis- 
sion originate from the explosion's O-burning layers. The regions where Channels 
2 (4.5 /mi) and 3 (5.6 /im) are strongest relative to Channel 4 show a spectrum 
that rises gradually to 21 /im, and then flattens or rises more slowly to longer 
wavelengths, along with higher ratios of [Ne II] to [Ar II]. We suggest that the 
dust and ionic emission in these locations arise primarily from the C- and Ne- 
burning layers. There are no bright lines in the Channel 3 spectra themselves, 
although strong line emission, perhaps from [Fe II], Br a or CO, must contribute 
significantly to the Channel 2 brightness. 

All of these findings are consistent with asymmetries deep in the explosion, 
producing variations in the velocity structure in different directions, but generally 
preserving the nucleosynthetic layering. At each location, the dust and ionic lines 
in the mid-infrared, and the hotter and more highly ionized optical and X-ray 
emission are then dominated by the layer currently encountering the reverse shock 
in that direction. 

Subject headings: supernova remnant: Cassiopeia A - dust - nucleosynthesis 

1. Introduction 

Cassiopeia A (Cas A) is the youngest supernova remnant (SNR) in our galaxy, thought 
to be the result of either a type lb or Iln supernova explosion (Chevalier & Oishi 2003) 
occurring in 1671 (Thorstensen, Fesen, & van den Bergh 2001) at a distance of 3.4 kpc 
(Reed et al. 1995). The remnant has been studied extensively at many wavelengths, and 
is one of the brightest radio and X-ray sources in the sky. Its primary structures are a 
105" radius bright ring surrounded by a 150" radius low surface brightness plateau (Braun 
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1987). The outer plateau is bordered by a thin X-ray ring identified as the outer shock in 
the circumstellar medium (CSM), with the broader, brighter interior ring originating from 
stellar ejecta that have encountered the reverse shock (Gotthelf et al. 2001). 

The X-ray emission is characterized by a thermal spectrum containing emission lines 
from highly ionized atoms. Optical emission from the remnant is dominated by chemically- 
enriched knots. Infrared emission was previously known to contain thermal continua from 
heated dust, line emission from ionized atoms, and, shortward of about 5 jum, synchrotron 
emission from electrons accelerated in shock regions (Jones et al. 2003; Rho et al. 2003). 
Submillimeter observations of Cas A prove difficult due to a molecular cloud complex along 
the line of sight and the presence of cold dust in the remnant is therefore still in question 
(Dunne et al. 2003; Krause et al. 2004; Wilson & Batrla 2005). The radio emission from Cas 
A is synchrotron radiation (Ginzburg & Syrovatskii 1965). 

Cas A's structure and dynamics reflect different progress into the Sedov- Taylor evolu- 
tionary phase in different directions (DeLaney & Rudnick 2003). The bright ring remains 
illuminated as new, successively slower-moving ejecta encounter the reverse shock and are 
heated and ionized (DeLaney et al. 2004). A Doppler analysis of the X-ray gas, studies of 
the optical knots, and the large abundance ratio of 44 Ti/ 56 Ni support an asymmetric ex- 
plosion (Reed et al. 1995; Hwang & Laming 2003; Willingale et al. 2003; Nagataki et al. 
1997). The MIPS images (Hines et al. 2004) showed both the main X-ray jet (Hwang, Holt, 
& Petre 2000) and a counterjet (Hwang et al. 2004), providing further evidence for explosion 
asymmetry. 

The progenitor of Cas A is generally believed to have been a WN star (i.e., a Wolf- 
Rayet star with high nitrogen abundance), due to the high abundances of N and H in some 
of the Fast-Moving Knots (FMKs, Kamper & van den Bergh 1976; Fesen & Becker 1991). 
The hydrodynamical model of Perez-Rendon, Garcia-Segura, & Langer (2002) suggests a 
29 — 30 M Q progenitor, while Young et al. (2006) find the overall data are best fit by 
a 15 — 25 M progenitor that loses its hydrogen envelope in a binary interaction. The 
pre-supernova wind produced a dense, clumpy medium (Chevalier & Oishi 2003) which is 
currently being shocked by the blast wave. The highest density shocked clumps are seen in 
optical emission as slow-moving Quasi-Stationary Flocculi (QSFs, van den Bergh & Kamper 
1985). 

Rich in optical, X-ray and infrared emission lines from ionized atoms (Fesen at al. 
2001; Douvion, Lagage, & Cesarksky 1999; Hwang, Holt, & Petre 2000), Cas A provides an 
important window into both quasi-equilibrium (pre-explosion) and explosive nucleosynthesis. 
Each layer contains several different elements, with C and O produced in He-burning, Ne 
and Mg first appearing through C-burning, and O and Al added with Ne burning (Woosley 
& Weaver 1995; Woosley & Janka 2005). When O and Mg burn, the heavier elements, Si, S, 
Ar, Ca are then produced and their burning products yield the Fe group elements. In Cas A, 
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much of this layered nucleosynthetic structure has been preserved following the supernova 
explosion, e.g., with layers of nitrogen- sulfur- and oxygen- rich ejecta seen beyond the outer 
shock (Fesen 2001). There is also some evidence for mixing of the nucleosynthetic layers 
on various scales. Optical "mixed emission knots" show both N and S lines, suggesting that 
high speed clumps of S ejecta penetrated through the outer N-rich layers (Fesen 2001). At 
infrared wavelengths, ISOCAM observations showed the presence of Ar and S in strong Ne 
knots, although the Ne and silicate emissions appeared anti-correlated (Douvion, Lagage, 
& Cesarksky 1999). The X-ray line data show some large-scale overturning of ejecta layers 
(Hughes et al. 2000; Hwang & Laming 2003); iron in the SE region is further out and is 
moving faster than the Si/O regions. 

The variations in ionic species at different locations, as probed by optical, infrared, and 
X-ray observations, are sensitive to temperature, density and ionization state. Thus, multiple 
wavelength detections of a species from different ionization states can help separate these 
effects from actual abundance variations (Vink et al. 2001). This provides one important 
motivation for our Spitzer studies. Density information can be derived from comparing 
different line strengths from the same ionization state, e.g. for [S III] (Houck et al. 1984). 
All of this information aids in reconstructing a picture of the inhomogeneities in the explosion, 
produced, e.g., through instabilities deep in the core (Foglizzo 2002; Blondin et al. 2003). 

Another key motivation for the Spitzer observations presented here was to understand 
the production and destruction of dust in Cas A. One possible location for dust production 
is in optically thick clumps of ejecta, such as proposed by Lucy et al. (1991) for SN1987A. 
They suggest that most of the dust is contained in such optically thick knots, with small 
grains being distributed diffusely between them. Ejecta knots that are dense enough should 
remain largely intact with the passage of the reverse shock and blast-wave. If the dust 
contained in dense clumps is not in equilibrium with the more diffuse X-ray gas, then it may 
remain at a colder temperature. 

In Cas A, Dwek et al. (1987) detected a strong excess of emission in Infrared Astronomy 
Satellite (IRAS) observations at 12 /iin to 100 ^m, largely from dust swept up by the 
supernova blast wave (the outer shock). Spitzer MIPS observations of Cas A (Hines et 
al. 2004) also detected this thermal dust emission from shocked circumstellar material. In 
addition, emission has been seen from a hot dust component associated with both the optical 
(Dwek et al. 1987; Fesen at al. 2001) and the X-ray ejecta (Douvion, Lagage, & Pantin 2001). 
Dust continua found in Cas A with the Infrared Space Observatory (ISO) were fit at 21 fim 
with proto-silicates (Arendt et al. 1999) or MgSi03 and Si02 (Douvion, Lagage, & Pantin 
2001). Those data suggest that the dust is heated continuously, presumably by the hot X-ray 
emitting gas (Dwek et al. 1987; Hines et al. 2004). As noted above, the amount of cold dust 
associated with Cas A is still uncertain. 

In this paper, we report the results of Spitzer Space Telescope images made using the In- 
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frared Array Camera (IRAC, Fazio et al. 2004) with brief supporting data from the Infrared 
Spectrograph (IRS, Houck et al. 2004) and images from other wavelengths. We show evi- 
dence for different nucleosynthetic layers currently encountering the reverse shock in different 
directions. In subsequent papers, we will address the physical conditions and dynamics of the 
gaseous material, and the detailed composition, temperature structure and mass estimates 
of the dust components. 

2. Observations 

IRAC observations covered the entire Cas A supernova remnant, including the outer 
shock, jet and counterjet regions. The IRAC images utilize four wide filters with central 
wavelengths of 3.6, 4.5, 5.6 and 8 fim for Channels 1, 2, 3 and 4, respectively. The data 
were taken on January 18, 2005. The observing strategy combined a mapping grid and 
dithers to yield a depth of coverage of at least 18 pointings over the entire remnant, with 
higher coverage in some overlap regions. At each pointing, a 0.6 and a 12 second frame was 
taken. The IRAC images have an angular resolution of 2-2.5" in Channels 1-2 and « 3" in 
Channels 3-4. The data were processed with the Sll version of the IRAC pipeline (Lowrance 
et al. 2006). The four IRAC images are shown in Figure I. 

The IRS Spectrograph was used on January 13, 2005 to spectrally map the full remnant 
(with portions of the outer structures missing from some slits), covering 5-15 /im (Short- 
Low module, SL) and 15-38 /im (Long-Low module, LL). Each included the two orders of 
wavelength, and the two long, low-resolution slits provided resolving powers of 64-128. The 
long- wavelength (15-38 [im) spectra were taken in a single large map with 4x91 pointings, 
using a single 6 second ramp at each position. To achieve the spatial coverage with the 
short-wavelength (5-15 fim) slit, a set of four quadrant maps were made, two with 4 x 87 
pointings and two with 3 x 87 pointings, using a 6 second ramp at each position. The 
mapped area ranged from 6.26' x 5.86' (SL) to 11.0' x 7.79' (LL), with offsets between 
the maps produced in each of the two orders in each module of 3.2' (LL) and 1.3' (SL), 
along the slit direction. The effective overlap coverage of all modules and orders is 4.9' by 
5.8'. The illustrative spectra presented here were processed with the S12 version of the IRS 
pipeline, using the Cubism package (Smith et al. 2006) to reconstruct the spectra at each 
position. They were extracted from the IRS cubes using areas from 10" to 33" across. No 
detailed matching of spectral overlap amplitudes between the SL and LL detectors was done. 
Occasional instrumental problems, such as spectral fringing at the short wavelength end of 
LL2, are not addressed and do not affect the analysis presented here. 

For comparisons with the IRAC images, we also performed near-infrared observations 
with a narrow Pa (3 filter, using the Palomar 200-inch Wide-infrared Camera (WIRC). The 
Pa (3 filter is cented on 1.282/im with a 1% width. The data were taken on Aug 15 and 16, 
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2005. The exposure time was nine times 90 seconds per sky position. 

3. Results 

3.1. IRAC images and comparisons to other bands 

The four IRAC images in Figure 1 each show the same overall structure of the remnant, 
including the bright ring, the surrounding low surface brightness plateau, and the eastern jet. 
The plateau region and the internal filamentary emission are most apparent in Channel 1. 
The large oval ring covering the northern third of the remnant is prominent in Channels 2 and 
4, weak in Channel 3, and just visible in Channel 1. Channels 3 and 4 show significant diffuse 
and patchy emission beyond the plateau, likely associated with the surrounding medium. 

A color image combining all four IRAC images is shown in Figure 2. There are large 
spatial variations in the relative strengths of the IRAC channels, resulting in the broad 
range of IRAC "colors". In order to examine the colors more quantitatively, we isolated 
those regions where each channel was strongest with respect to Channel 4, and determined 
the mean surface brightness in those regions for each of the four channels. The results are 
shown in Figure 3. The very large jump to Channel 4 (ps 10) is likely caused by the presence 
of [Ar II] and [Ar III] emission in the IRAC band, as discussed further below. Unfortunately, 
the IRS spectra do not cover Channels 1 and 2, and cover only part of Channel 3, so we 
cannot perform a quantitative analysis of these various IRAC colors; below, we suggest a 
few possible contributors to Channels 2 and 3. 

In order to understand the origins of these IRAC color differences, we therefore first 
compare the images in each IRAC channel with those from other bands. The IRAC Channel 
1 image is very similar to that seen in the radio (Figure 4), with interior, ring, and plateau 
emission. By contrast, the bright ring dominates the emission in Channel 4, mirroring 
Cas A's appearance in the optical, in the 24 /xm MIPS image, and in X-ray line images 
(Figure 5). Channels 2 and 3 also are dominated by bright ring emission, although there 
are distinct differences in the brightness of various features between them and Channel 4. 
Figure 6 shows a comparison of Channels 2 and 3 with the [Fe II] (18 /im) and Pa (3 (1.3 //m, 
with some [Fe II] contamination) images. The results of these comparisons are all discussed 
in more detail below, after we have examined the correspondence between the IRAC colors 
and the shapes of the IRS spectra. 

We find no evidence for emission in the IRAC images from Cas A's compact X-ray 
source (Tannanbaum 1999) against the variations in flux near the center of the remnant. 
The 3-sigma upper limits are 50 ^Jy for Channels 1 and 2, 100 /iJy for Channel 3, and 
220 ^Jy for Channel 4. 
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3.2. IRS Spectra 



The 5.3-38.5 /im IRS spectrum for the full remnant is presented in Figure 7, showing 
both the average continuum shape as well as Doppler-broadened ionic line emission. We 
indicate where the IRS spectral coverage overlaps IRAC Channels 3 and 4; there is no IRS 
coverage in the IRAC Channels 1 and 2 bands. We also extracted spectra from 22 different 
regions in the remnant, chosen to explore a broad range of possible physical properties by 
comparing the IRAC color image to MIPS 24 /jm, X-ray, optical, and radio images. These 
sample spectra showed a variety of relative line strengths and continuum shapes, especially in 
the relative strength of the peak around 21 //m, as expected from the ISO work of Douvion, 
Lagage, & Pantin (2001). We found that the spectra fell into three major categories, as 
follows: Broad - showing a gentle peak around 10 //m and rising to a very broad, gradual 
peak around 26 /im, with little line emission; Strong 21 [im - showing a 2-3 /im wide strong 
asymmetric peak at 21 /im, similar to those studied with ISO (Douvion, Lagage, & Pantin 
2001), along with strong lines of Ar, Ne, Si, S, and 26 /im Doppler-blended Fe and O; and 
Weak 21 fim - rising gently through 21 /im and gradually becoming shallower to longer 
wavelengths, accompanied by stronger [Ne II], but relatively weaker [Ar II] lines. Figure 7 
shows the average shape for each of these classes; the variations of shape within each class 
are indicated by the grey bands. 

We find a good correspondence between the IRS spectral shapes and the IRAC colors, 
as illustrated in Figure 8. When Channel 4 (with major contributions from [Ar II]) is strong 
with respect to Channel 2, the IRS spectra show the strong 21 /im shape. When Channel 
4 is weaker, the ratio of Channels 2 and 3 distinguishes between broad spectra and weak 
21 /im spectra. A relatively high ratio of Channel 1 to Channels 2 or 3 is also a good indicator 
of broad spectra. Using the various IRAC colors, we can identify locations in the remnant 
where each IRAC channel in turn appears strongest compared to the other channels. Table 
1 summarizes the typical properties of regions where each respective channel is so enhanced. 
We now look at each of the IRAC channels in turn, and discuss the likely origins of their 
emission. 



3.3. Strong IRAC Channel 1 (3.2 - 3.9 /im) 

The Channel 1 image is shown in Figure 4, along with a A 6cm radio image from DeLaney 
(2004). The detailed correspondence between these two images shows that Channel 1 is 
dominated by synchrotron emission. The synchrotron nature of the emission at 2.2 /im was 
first suggested by Gerardy & Fesen (2001) based on its morphological similarity to the radio 
emission, and then established both by its polarization (Jones et al. 2003) and brightness 
at levels expected from extrapolations of the radio spectrum (Rho et al. 2003; Jones et al. 
2003). We now extend the detection of synchrotron radiation to the mid-infrared; in IRAC 
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Channel 1 we see the same bright ring, faint plateau and filamentary structures as in the 
radio, at brightness levels comparable to those calculated from an extrapolation of the radio 
spectrum. In the forward shock region, where IRAC Channel 1 is most enhanced relative to 
the other channels, synchrotron radiation makes substantial contributions to all the IRAC 
Channels (Figure 3). This region also shows significant 4-6 keV X-ray emission in the form 
of a thin rim of emission at the edge of the radio plateau. The X-ray rim has been identified 
as marking the location of the outer (forward) shock (Gotthelf et al. 2001) and is likely 
dominated by synchrotron emission (Vink et al. 1999; DeLaney et al. 2004). A detailed 
analysis of the radio/infrared/X-ray spectral shape holds important clues to the relativistic 
particle acceleration mechanism (Ellison, Decourchelle & Ballet 2005) but is beyond the 
scope of this paper. 

Substantial emission from the forward shock region can also be seen in the MIPS 
24 fim image (Figure 5, Hines et al. 2004). However, the 24 /im brightness reaches up 
to a factor of ^50 above the values extrapolated from IRAC Channel 1, so it must be due 
to another emitting component. IRS spectra of this spatial component show the charac- 
teristic "broad" shape shown in Figure 7. Previous observations with ISOCAM (Arendt et 
al. 1999) found these spectra rising to 18 /zm; with IRS, we now see a broad peak around 
26 //m, a smaller "bump" around 9 /im, and little or no line emission. The IRS spectrum 
of the forward shock region shown here can be approximately fit by a blackbody Planck 
function with a temperature of 113 K, multiplied by the absorption efficiency calculated 
for grain models for R v = 3.1 (Weingartner & Draine 2001). The silicate emission feature 
between 9-11/im from the interstellar medium can also be seen; such silicate dust from the 
circumstellar /interstellar medium is expected in the forward shock region. 



3.4. Strong Channel 4 (6.5 - 9.4 /mi) 

The IRAC Channel 4 image (Figures 1 and 5 ) contains both continuum and line emis- 
sion, with [Ar II] at 6.99 /im contributing approximately 40% of the brightness, while [Ar III] 
at 8.99 /im contributes «3%. The 5-38 /im IRS spectra of Channel 4 dominated regions 
also show significant emission from [Ne II] (12.8 /im) and [S IV] (10.5 /im). Strong lines 
from either [Fe II] and/or [O IV] (26 /im), [Si II] (34.8 /im) and [S III] (33.5 /im), are also 
seen, but these are not exclusive to bright Channel 4 regions. In most locations associated 
with the bright ring, Channels 2 and 3 are weak with respect to Channel 4 (e.g., as seen in 
Figure 3). 

The shape of the continuum of bright Channel 4 regions is characterized by a 2-3 /zm 
wide peak at 21 /xm, similar to that seen by Arendt et al. (1999) who modeled it as due 
to magnesium proto-silicates. This 21 /im peak thus leads to an excellent correspondence 
between Channel 4 and the MIPS 24 /im image (Figure 5). However, the MIPS image also 
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shows strong emission from the outer shock, whose spectra peak around 26 ^m, without any 
corresponding strong Channel 4 emission. 

Many of the details of the Channel 4 emission can also be seen in the optical HST 
WFPC2 images; Figure 5 shows [S II], [O II] and [O I] emission from dense shocked ejecta 
in the F675W image (Fesen at al. 2001). Some of the large-scale differences between the 
Channel 4 and WFPC2 images are due to variations in optical extinction. The 0.3 - 10 keV 
X-ray emission (Hwang, Holt, & Petre 2000) also shows some similarities to the Channel 4 
image; a better correspondence is seen in the X-ray Si (shown in Figure 5) and S emission. 
Ionic lines from lower ionization states of these two elements are seen in the IRS spectra of 
strong Channel 4 regions. 

The ratio of [Ar II] to [Ne II] is highest in the bright Channel 4 regions, partly due to 
selection. In order to quantify this, we created continuum-subtracted line maps around the 
[Ne II] (12.8 //m) and [Ar II] (6.99 /itm) images, and calculated the average ratio of those 
lines (see Table 1) in regions where the Channel 4 brightness was above 25 MJy/sr. This 
will serve as a standard of comparison for the [Ar II] /[Ne II] ratio in places where the other 
IRAC channels are strongest relative to Channel 4. 

3.5. Relatively Strong IRAC Channels 2 (4 - 5 /im) and 3 (5 - 6.4 fim) 

Channel 2, and to a lesser extent, Channel 3, have significant contributions from syn- 
chrotron radiation in some locations such as the northern forward shock region and interior 
filamentary structures (see Figure 3). We have therefore subtracted the Channel 1 image 
from Channels 2 and 3 to create the residual images in Figure 6. The similarities between 
Channel 2 and Channel 4 are now more apparent, although the relative brightnesses of fea- 
tures vary by at least a factor of 5. Channels 2, 3, and 4 all trace out the same large oval 
structure in the North, e.g., but their relative strengths vary abruptly between nearby knotty 
structures (see Figure 2). Sometimes, as in the far north of the bright ring, these changes 
are actually due to dynamically distinct but superposed features. 

Bright Channel 2 regions show both strong 21 fj,m and weak 21 pm spectra, but when 
Channel 2 is brightest with respect to Channel 4, we find only weak 21 /mi spectra, along 
with low [Ar II] /[Ne II] ratios (see Table 1). These low ratios occur because [Ar II] (Channel 
4) is weak in these regions. Most of the relatively bright Channel 3 regions also have low 
[Ar II] /[Ne II] ratios and weak 21 fj,m spectra, although some show modest peaks around 
21 /xm. 

There is no IRS coverage for Channel 2; we consider possible line contributions below. 
Channel 3 is covered by the IRS from 5.3 to 6.4 /im, but there is no coverage between 5.0 
and 5.3 jum. In the accessible wavelength range we found no ionic lines dominating the 
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Channel 3 brightness, although there are very weak contributions in some locations from 
[Fe II] 5.3 /tm at the edge of the IRS band. To see where [Fe II] might be important for 
Channel 3, we present the 18 /tm [Fe II] image in Figure 6. [Fe II] emission is also seen 
at 1.64 /mi (Rho et al. 2003) and in the spectra of some FMKs around 1.2 (Gerardy & 
Fesen 2001). Comparison with our near infrared measurements suggests that [Fe II] might 
contribute up to 50% of the Channel 3 emission in some isolated locations, but a negligible 
amount in many bright Channel 3 regions. The rest of Channel 3 emission is likely from the 
dust continuum. When the synchrotron emission is subtracted from Channel 3, we also find 
a few isolated bright patches in the northeast jet and elsewhere (Figure 6). These bright 
patches are not coincident with features seen in Channels 2 and 4, and their IRS spectra 
show the "broad" shape characteristic of the forward shock. 



3.5.1. What is Channel 2 (4 - 5 fim) ? 

In the absence of IRS coverage in the Channel 2 wavelength range, we summed the 17 
available SWS spectra from the ISO archives (http://www.iso.esac.esa.int/ida) but found 
no strong lines shortward of the 6.99 /tm [Ar II] line. However, the fact that Channel 2 is 
usually brighter than both Channel 1 and Channel 3 (Figure 3) indicates the presence of 
line emission. 

There are two distinct questions regarding the origins of the Channel 2 emission - what 
dominates when Channel 4 (and [Ar II]) is also strong, and what dominates when Channel 
4 (and [Ar II]) is weaker, but [Ne II] is still strong? [Fe II] has several lines in the Channel 
2, 4 - 5 /mi, band, but the 5.3 /mi emission is quite weak. Looking at the 18 /tm [Fe II] 
structure, there is little or no emission where Channel 2 is strong relative to Channel 4 (see 
Figure 9 for the locations of these regions), so [Fe II] is unlikely to provide the missing lines 
in Channel 2. 

If pieces of the hydrogen envelope survived the WR-wind stage of Cas A, then Br a at 
4.05 /tm could be present in the ejecta. Since ground-based spectroscopy of Br a is difficult, 
we obtained an image of the Pa (3 1.28 /im line using the Wide-Field Infrared Camera on 
the 5m Hale Telescope at Palomar Observatory. Some nearby [Fe II] lines, seen in FMKs 
by Gerardy & Fesen (2001) also fall into this filter. In a number of regions, the Pa (3, [Fe 
II] and Channel 2 images trace out the same structures, so Br a might be responsible for 
some Channel 2 emission. These are also the regions where Channel 4 is strong. However, 
there are other locations where Channel 2 is strong, such as the jet and the crescent shaped 
regions seen in Figure 9, that are weak or absent in the Pa (3 image, so Br a is unlikely to 
be playing a key role. We also see no evidence for Pfa at 7.46 /tm. 

Another strong candidate is the CO fundamental bandhead around 4.76 /mi, such as 
seen in regions with shocked CO (Gonzlez- Alfonso et al. 2002). This bandhead has been 
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detected in SN 1987A (Meikle et al. 1989; Kotak et al. 2005), as well as the first overtone 
at 2.29 fxm (e.g., Catchpole et al. 1987), so we know that CO can form in supernova ejecta. 
The resolution of this issue requires sensitive spectra in the 2-5 //m band. 

Another possibility is H 2 , which we see in our IRS images around 17 ^m. However, the 
H 2 is largely exterior to the remnant, possibly associated with the surrounding CO clouds 
(Liszt & Lucas 1999). If H 2 were dominant in Channel 2, it should also be strong both in 
5-6.4 fim and 6.4-9.4 fxm spectra, but is not seen. 

Hell 8-7 line recombination occurs at 4.76 fim. If this were responsible for Channel 2 
emission, we should also have Hell 9-8 emission at 6.95 fim and 10-9 emission at 9.71 /j,m. 
The former is unfortunately coincident with the extremely bright [Ar II] emission, and we 
find no evidence for the latter anywhere in the remnant. At present therefore, the origin of 
the line emission in Channel 2 is unclear. 



4. Discussion 

We have found the mid-infrared radiation from Cas A to arise from a number of different 
components - at short wavelengths, synchrotron radiation and at longer wavelengths, low 
ionization lines from Ne, O, Si, S, Ar and Fe ejecta, and shock-heated dust from both ejecta 
and CSM. The ejecta at different locations are further distinguished from each other by their 
colors in the IRAC bands, by the relative line strengths of different elements, and by the 
shape of their dust continua. The variations occur on both small and large spatial scales. The 
same spatial variations characterize the optical and X-ray line emission from ejecta, although 
these are from much higher ionization states. In this discussion, we briefly summarize the 
structure of the multiwavelength appearance of the ejecta, and the implications of these new 
Spitzer observations for the dynamics of Cas A's explosion. 

A consistent picture of the ejecta structure emerges from images at different wavelengths. 
Line emission from elements such as Si and S dominate the optical (Fesen at al. 2001) and X- 
ray (Hwang, Holt, & Petre 2000) emission, from moderate and high ionization states. These 
appear structurally as a partially illuminated bright circular ring at the same location as the 
radio bright ring. In addition, both optical and X-ray observations show the NE "jet" and 
two interior elliptical rings towards the North. IRAC Channel 4 (6.4 - 9.4 yum), which has 
major contributions from [Ar II] and [Ar III] emission, shows all of these same structures. 

A quite different picture of the ejecta emerges from IRAC Channel 3 (5 - 6.4 /im), which 
is dominated by dust, and in regions where Channel 2 (4 - 5 ^m) is strongest with respect 
to Channel 4. Regions of high (Channel 2 / Channel 4) can be seen in Figure 9, where 
the most prominent features are two bright crescent features. The ratio of [Ar II] /[Ne II] 
from our IRS spectra (see Table 1) is low in these regions, so they are relatively Ne-rich, 
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although Ne does not directly contribute in the IRAC bands. The northern crescent is also 
shown in Figure 9 overlaid on the WFPC2 F450W image, showing that the same crescent 
structure appears in [O III] AA 4959, 5007 emission. This feature is also seen in an X-ray 
image in the oxygen emission around 0.64-0.71 keV, using unpublished ACIS data from our 
proper motion studies (DeLaney et al. 2004). The brightness of the northern crescent in the 
F850LP WFPC2 image, which is primarily sensitive to [S III] AA 9069, 9531 emission, is 
a factor of two lower, compared to the [O III] emission, than in surrounding regions. The 
southern crescent falls in a very distinct X-ray gap, as seen in Figure 9. Multiple epoch HST 
images show brightening [O III] emission at this location (R. Fesen, private communication). 
Although there is some X-ray Si, S, or Fe emission in the general area of the crescents 
(Hwang, Holt, & Petre 2000), they do not show the detailed structural correspondence seen 
in the O X-rays. 

The ejecta thus appear to spatially segregate in a way consistent with different nucle- 
osynthetic layers. The brightest optical, X-ray and 8 fim (Channel 4) emission is dominated 
by lines from the O-burning layers, e.g., Si, S and Ar. Other regions, such as the crescents, 
are distinguished by their relatively strong 4-6 /im emission and the relatively strong 
products of C-burning, neon (infrared) and oxygen (optical and X-ray). This suggests that 
the appearance of different elements at different locations in the remnant reflects which 
nucleosynthetic layer is locally illuminated. 

This picture receives strong support from the different types of dust found associated 
with regions of different IRAC colors. Where Channel 4 is strong, the IRS spectra show a 
broad spectrum with a 21 /im peak, which requires the presence of silicates, e.g. MgSiOs 
or Fe or Mg proto-silicates (Lagage et al. 1996; Arendt et al. 1999; Rho et al. 2006). Thus, 
the O-burning product Si is also seen in the dust. Where Channels 2 and 3 are strongest 
relative to Channel 4, the gently rising weak 21 /iin spectra are dominated by A1 2 3 , i.e., the 
C-burning products. The detailed models of dust temperature and composition, including 
contributions from carbon dust, are presented in Rho et al. (2006). 

These findings provide a new perspective on the distribution of different elements on 
Cas A's sparsely-covered spherical shell, which appears in projection as the bright ring. Ap- 
parent differences in composition at different locations could have appeared from variations 
in temperature and/or ionization states. Similarly, variations in ionization timescale (Hwang 
& Laming 2003) and density can significantly affect which elements are seen. However, the 
results presented here indicate that differences in apparent composition likely reflect the 
actual local composition. In some locations, we find multiple indicators for only C-burning 
products, e.g. neon, and oxygen in the form of [O III], [O VIII] and AI2O3 dust. In other 
locations we see O-burning products. For example, when sulfur is seen, e.g., it appears as 
[S III] (33 fim, see Fig. 7) and He- and H- like S (X-ray S XV and S XVI lines between 
?a2.4 and 3.1 keV). Silicon, where it is present, appears as strong 21 /im dust (i.e., silicates), 
[Si II] (35 H. and He ~ and H- like Si ( x - r ay Si XIII and Si XIV lines between ml. 8 and 
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2.6 keV). These variations in composition at different locations likely reflect asymmetries in 
the original explosion. 

We briefly outline a simple dynamical model for Cas A which incorporates these findings. 
The blast wave from the explosion was quite symmetric, as seen by the nearly circular 
appearance of the outer and global reverse shocks (Gotthelf et al. 2001), and has now swept 
up sufficient mass to decelerate the outer shock by a factor of ssl.5 (DeLaney et al. 2004) 
from free expansion. This produces an inward-moving (in the frame of the explosion) global 
reverse shock in the diffuse X-ray gas. Clumps of ejecta traveling at ^5000 km/s will be 
encountering this reverse shock at the current epoch, driving a local reverse shock back into 
the clumps. This heats and compresses them, making them visible optically (Morse et al. 
2004). Stripping, heating, ionization and disruption of these clumps then lead to decelerated 
X-ray and radio-emitting features (DeLaney et al. 2004; Anderson et al. 1994) and eventual 
disappearance. The reverse shock then is a slowly moving front that is successively overtaken 
and illuminated as a bright ring, by increasingly slower-moving undecelerated ejecta from 
the initial explosion. An illustration of the key features of this picture is shown in Figure 10. 

Superposed on this symmetric structure are major variations, the most prominent being 
the jet and counterjet regions. They arise deep in the explosion, producing fast-moving S- 
rich ejecta (Fesen et al. 2006), as well as emission from Si group elements (Hwang et al. 
2004). Roughly perpendicular to this axis, we now find evidence for a much slower moving 
bipolar structure, in the form of the crescents of C-burning material in the infrared, optical, 
and X-ray. In order to reach the global reverse shock and become visible now, this material 
must have moved at a free expansion speed of ~5000 km/s. However, along this axis only 
these upper layers are now encountering the reverse shock. Material from the O-burning 
layers is not seen either at the reverse shock, or in the fast-moving outlying knots in these 
directions (Fesen 2001). 

This scenario leads to the suggestion that if we could wait sufficiently long, we would 
see the O-burning layers encounter the reverse shock along the crescents' axis, and become 
visible at all wavelengths. We find evidence that this actually may be occurring, because of 
the presence of [Si II] and [S III] emission interior to and somewhat overlapping the bright 
ring. 

5. Conclusions 

1. The four IRAC bands dominate in different regions of the Cas A supernova remnant, 
echoing structures seen in optical, X-ray and radio images. 

2. IRAC Channel 1 is dominated by infrared synchrotron radiation; where Channel 
1 dominates, the broadband spectra have a distinct shape, gently peaking around 26 fim, 
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which we attribute to forward-shock- heated circumstellar dust. 

3. IRAC Channel 4 has a significant contribution from both [Ar II] emission and con- 
tinuum. Where Channel 4 dominates, the dust continuum peaks strongly around 21 /an, 
signifying the presence of silicates. 

4. Where IRAC Channels 2 and 3 are strongest with respect to Channel 4, [Ar II] is 
weaker relative to [Ne II]. The continuum in these regions rises slowly or levels off at 21 /im, 
showing the absence or only weak presence of silicate dust. 

5. The relatively strong Channel 2,3 regions show optical and X-ray oxygen emission 
and an absence of silicon and sulfur. 

6. We present a dynamical picture for the Cas A's explosion in which cones of ejecta are 
produced with different velocities in different directions, In some directions, only the upper 
C-burning layers have reached the reverse shock, while in other directions, deeper O- and 
Si- burning layers have done so. A particularly interesting bipolar structure, perpendicular 
to the jet/counterjet axis, has been identified. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided by NASA through an award issued 
by JPL/Caltech. 
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Fig. 1. — IRAC images of Cassiopeia A. Channel 1: 3.2 - 3.9 /im. Channel 2: 4.0 
5.0 yum. Channel 3: 5 - 6.4 yum. Channel 4: 6.5 - 9 /im. The bright ring 
marks the location of the large-scale reverse shock, and the outer edge of the faint, dif- 
fuse plateau is the location of the forward shock. For full-resolution images, please see 
http:/ /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Fig. 2. — Color combination of IRAC images. Channel 1: blue; Channel 2: green; 
Channel 3: orange-yellow; Channel 4: red. For full-resolution images, please see 
http: / /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Fig. 3. — Broadband IRAC spectra of locations in Cas A where each of the respective 
IRAC Channels are enhanced relative to the other channels. An estimate of the synchrotron 
contribution is shown as a grey line superposed on the Channel 1 enhanced spectra. For 
the other spectra, only Channel 1 can have a dominant synchrotron component because the 
brightness rises much too rapidly at longer wavelengths. For full-resolution images, please 
see http:/ /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Fig. 4. — Left: IRAC Channel 1; Right: VLA radio image at 6cm from 
DeLaney (2004), 1.3" resolution. For full- resolution images, please see 

http:/ /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Fig. 5. — Comparisons of IRAC Channel 4 with images from other wavelengths - 
MIPS 24 jum (Hines et al. 2004), HST WFPC2 F675W, (Fesen at al. 2001), and 
Chandra ACIS X-rays for the Si line around 1.85 keV, using the epoch 2002 proper 
motion data set of DeLaney et al. (2004). For full-resolution images, please see 
http:/ /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Fig. 6. — Synchrotron subtracted images of IRAC Channels 2 (top) and 3 (bottom) for 
comparison to possible indicators of lines in 4 - 5 ^m range - Pa (3, (with some [Fe II] 
contamination) at wl.3 /mi, and [Fe II] image at ~18 /tm. For full- resolution images, please 
see http:/ /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Wavelength (microns) Wavelength (microns) 

Fig. 7. — IRS spectra of the total remnant and the three major spectral classes. The spectra 
show the mean of the class, and the grey regions show the rms scatter within the class. Within 
each class, all spectra are normalized to have the same brightness at 30 /im. Small (up to ~ 
10%) cosmetic adjustments have been made to normalize the SL and LL brightness scales. For 
full- resolution images, please see http://webusers.astro.umn.edu/~jennis/iracpaper.html. 
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Fig. 8. — IRAC colors for the three major classes of IRS spectra. For full-resolution images, 
please see http:/ /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Fig. 9. — Channel 2 (red) and Channel 4 (green). The two red crescent regions (rel- 
atively high [Ne II] /[At II]) are shown on right, overlaid on images from other wave- 
bands. Top: Contours of Ch2/Ch4 ratio overlaid on the WFPC F450W image ([O III] 
AA4959, 5007) Fesen at al. (2001). Bottom: Contours of Ch2/Ch4 ratio overlaid on 
Chandra Si image (Hwang, Holt, & Petre 2000). For full-resolution images, please see 
http: / /webusers. astro. umn.edu/~jennis/iracpaper. html. 
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Fig. 10. — Illustration of dynamical model described in text. In different directions, different 
layers from the explosion are currently reaching the reverse shock. There, they develop 
internal shocks and become visible for a short time across all wavebands before fading. 
X-ray and radio emission do not appear until significant deceleration has occurred. For 
full- resolution images, please see http://webusers.astro.umn.edu/~jennis/iracpaper.html. 



27 



Table 1: Characteristics of emission for locations where each IRAC channel is enhanced 
relative to the other channels. 



Channel 


Emission 


[Ar II] /[Ne II] 


Broadband shape 


Origin 


AA fxm 






(Fig. 2 color) 




1 


synchrotron 


no lines 


Broad 


Forward shocked 


3.2 - 3.9 






(blue) 


CSM 


2 


CO?, Br a? 


1.8 


Weak 21 /mi 


C-burning products 


4.0 - 5.0 


[Fe II]?, dust 




(green) 


(0, Si- burning 
products) 


3 


dust, [Fe II] 


1.7 


Weak 21 /mi 


C-burning products 


5.0 - 6.4 






(orange/yellow) 


(0, Si- burning 
products) 


4 


[Ar II], [Ar III] 


3.7 


Strong 21 /mi 


O, Si burning 


6.5 - 9.4 


dust 




(red) 


products 



Notes: The primary contributor to each IRAC channel is shown in boldface. 
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